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ABSTRACT 

We report the gravitational microlensing discovery of a sub-Saturn mass 
planet, MOA-2009-BLG-319Lb, orbiting a K or M-dwarf star in the inner Galac- 
tic disk or Galactic bulge. The high cadence observations of the MOA-II survey 
discovered this microlensing event and enabled its identification as a high mag- 
nification event approximately 24 hours prior to peak magnification. As a result, 
the planetary signal at the peak of this light curve was observed by 20 different 
telescopes, which is the largest number of telescopes to contribute to a planetary 
discovery to date. The microlensing model for this event indicates a planet-star 
mass ratio of q = (3.95±0.02) x 10~ 4 and a separation of d = 0.97537±0.00007 in 
units of the Einstein radius. A Bayesian analysis based on the measured Einstein 
radius crossing time, t-^, and angular Einstein radius, 9-^, along with a standard 
Galactic model indicates a host star mass of M L = 0.38tQ;fg M & and a planet 
mass of M p = 50±^ M ffi , which is half the mass of Saturn. This analysis also 
yields a planet-star three-dimensional separation of a = 2.4l ;g AU and a dis- 
tance to the planetary system of = 6.llJ' 2 kpc. This separation is ~ 2 times 
the distance of the snow line, a separation similar to most of the other planets 
discovered by microlensing. 

Subject headings: gravitational lensing: micro - planetary systems 

1. Introduction 

We present the eleventh microlensing planet, following ten previous discoveries (Bond 
et al. 2004; Udalski et al. 2005; Beaulieu et al. 2006; Gould et al. 2006; Gaudi et al. 2008; 
Bennett et al. 2008; Dong et al. 2009; Sumi et al. 2010; Janczak et al. 2010). Microlensing 
is unique among exoplanet detection methods in that it is sensitive to planets with masses 
down to 1 M ffi (Bennett & Rhie 1996) at relatively large separations, typically between 
1 AU and 6 AU, depending on the mass of the host star. These separations are generally 
beyond the "snow line" at ~ 2.7 AU M/M© (Ida & Lin 2004; Lecar et al. 2006; Kennedy 
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& Kenyon 2008), the region where planets can form most quickly, according to the core 
accretion theory. Microlensing confirms this expectation, as a statistical analysis of the 
prevalence of planets found by microlensing shows that Saturn-mass planets beyond the snow 
line are more common than the higher mass gas giants found by radial velocities in shorter 
period orbits (Gould et al. 2010), although the microlensing results are consistent with an 
extrapolation of the radial velocity results for solar-mass stars to larger orbital distances 
(Cumming et al. 2008). Furthermore, Sumi et al. (2010) have shown that the number of 
planets (per logarithmic interval) increases with decreasing mass ratio, q, as q-°- 7±0 - 2 ^ down 
to ~ 10 M®. So, cold Neptunes seem to be even more common than cold Saturns. While 
the number of planets found by microlensing is relatively small, it is the cold-Neptunes and 
Saturns discovered by microlensing that represent the most common types of exoplanet yet 
to be discovered. Microlensing has also found the first Jupiter /Saturn analog planetary 
system (Gaudi et al. 2008; Bennett et al. 2010), and it should soon be possible to use the 
microlensing results to determine how the properties of exoplanet systems vary with distance 
from the Galactic center. 

Searches for exoplanets via the microlensing method are currently conducted by two 
survey groups, the Microlensing Observations in Astrophysics (MOA; Bond et al. 2001; 
Sumi et al. 2003) and the Optical Gravitational Lensing Experiment (OGLE; Udalski 2003), 
which monitor ~ 40 deg 2 of the Galactic bulge to identify stellar microlensing events that 
can be searched for planetary signals. The planetary signals have durations that range from 
a few hours to a few days, so a global network of telescopes is needed to search for and 
characterize planetary signals. The follow-up groups that complete this telescope network 
are the Microlensing Follow- Up Network (/iFUN), RoboNet, Microlensing Network for the 
Detection of Small Terrestrial Exoplanets (MiNDSTEp), and the Probing Lensing Anomalies 
NETwork (PLANET). These narrow field-of-view follow-up telescopes can provide very high 
cadence observations of a small number of events that are known to be interesting, due to 
known or suspected planetary deviations in progress (Sumi et al. 2010) or high magnification 
events, which have very high planet detection efficiency (Griest & Safizadeh 1998; Rhie et 
al. 2000; Rattenbury et al. 2002). The very wide (2.2 deg 2 ) field-of-view of the MOA-II 1.8m 
telescope with 80M pixel CCD camera MOA-cam3 (Sako et al. 2008) provides high cadence 
survey observations of the entire Galactic bulge, and this allows MOA to identify suspected 
planetary deviations in progress and to predict high magnification (A max ^ 100) for events 
with relatively short timescales (Einstein radius crossing time < 20 days.). MOA-2009- 
BLG-319 is one such short timescale high magnification event that was identified as a high 
magnification event based on MOA data taken ~ 24 hours prior to peak magnification. 

In this paper, we report the discovery of a sub-Saturn mass planet in the microlensing 
event, MOA-2009-BLG-319. We describe the observations and data sets in Section 2. The 
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light curve modeling is presented in Section 3. We discuss the measurement of the source 
magnitude and color in Section 4, and derive the angular Einstein radius in Section 5. In 
Section 6, we search for a microlensing parallax signal. In Section 7, we use a Bayesian 
analysis to estimate the masses and distances of the host star and the planet, based on the 
angular Einstein radius and microlens parallax. We present our conclusion in Section 8. 

2. Observations 

For the bulk of the 2009 observing season, the MOA group was the only microlensing 
survey group in operation because the OGLE group completed the OGLE-III survey on 3 May 
2009, in order to upgrade to the OGLE-IV camera with a much wider field-of-view. Prompted 
in part by this fact, MOA adopted a new observing strategy for the 2009 observing season in 
order to increase the planet detection efficiency The top 6 fields (a total of 13.2 deg 2 ) yielded 
54% of the microlensing events found by MOA in previous seasons, and these were observed 
every 15 minutes. The next 6 fields (with 25% of the previous years' events) were observed 
every 47 minutes, and most of the remaining 10 fields were observed every 95 minutes. This 
new observing strategy yielded 563 microlensing alert events in 2009, an increase of about 
100 over the 2008 total. MOA-2009-BLG-319 was the first of four of these events to yield an 
apparent planetary signal. 

The event MOA-2009-BLG-319 [(R.A., decl.) J20 oo.o=(18' l 06 m 58M3, -26°49T0" .89), (/, 
6) = (4.202, -3.014)] was detected and announced as a normal microlensing alert event by 
the MOA collaboration on 20 June 2009 (HJD' = HJD - 2450000 = 5003.056). The dis- 
covery announcement provided a model for this event, which indicated that this was a 
high- magnification event, and so MOA immediately began follow-up observations in the / 
and ^-bands with the University of Canterbury's 0.6m Boiler & Chivens (B&C) telescope 
at Mt. John Observatory. Public access to the MOA photometry over the next two nights, 
led the /iFUN, RoboNet, and MiNDSTEp collaborations to begin observations of this event 
~ 2.5 days after its discovery. Three days after the discovery, the MOA data indicated that 
this event was quite likely to reach high magnification, and the /xFUN group issued a high- 
magnification alert by email to all interested observers, estimating a peak magnification of 
A max > 100 (at 1-cr ) 18 hours later at HJD' = 5006.875. This alert message noted "low-level 
systematics" in the MOA data, which were, in fact, not systematic errors at all. Instead, 
this light curve feature was the first (weak) planetary caustic crossing. Then, 14 hours 
later at June 24 UT 01:12 HJD' ~ 5006.55), data from the /iFUN SMARTS CTIO tele- 
scope in Chile provided clear evidence for a second, much stronger, caustic crossing feature, 
which prompted /iFUN to issue an anomaly alert. This feature was readily confirmed by the 
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MiNDSTEp observer at La Silla from data already in hand (see Fig. 1). A large number of 
telescopes responded to this anomaly alert, resulting in continuous photometric monitoring 
of the remainder of the planetary signal with no gaps larger than 5 minutes until after the 
planetary signal finished, some ~ 20 hours later. 

The complete data set for MOA-2009-BLG-319 consists of observations from 20 different 
observatories representing the MOA, /iFUN, RoboNet, MiNDSTEp, PLANET groups, as 
well as the InfraRed Survey Facility (IRSF) telescope in South Africa. Specifically, the data 
set includes data from the following telescopes and passbands: MOA-II (New Zealand) 1.8m 
wide-i?-band, the Mt. John Observatory B&C (New Zealand) 0.61m / and V bands, /iFUN 
Auckland Observatory (New Zealand) 0.4m i?-band, /iFUN Bronberg (South Africa) 0.35m 
unfiltered, /iFUN SMARTS CTIO (Chile) 1.3m V, I, and H bands, /iFUN Campo Catino 
Austral (CAO, Chile) 0.5m unfiltered, /iFUN Farm Cove (New Zealand) 0.35m unfiltered, 
/iFUN IAC80 (Tenerife, Spain) 0.8m I band, /iFUN Mt.Lemmon (Arizona, U.S.A.) 1.0m 
/ band, /iFUN Southern Stars Observatory (SSO, Tahiti) 0.28m unfiltered, /iFUN Vintage 
Lane Observatory (New Zealand) 0.41m unfiltered, /iFUN Wise (Israel) 0.46m unfiltered, 
/iFUN Palomar (U.S. A) 1.5m I band, RoboNet Faulkes Telescope North (FTN, Hawaii) 
2.0m SDSS-7 band, RoboNet Faulkes Telescope South (FTS, Australia) 2.0m SDSS-J band, 
RoboNet Liverpool Telescope (La Palma) 2.0m SDSS-J band, MiNDSTEp Danish (La Silla) 
1.54m I band, PLANET Canopus (Australia) 1.0m / band, PLANET SAAO (South Africa) 
1.0m I band, and IRSF (South Africa) 1.4m J, H and K§ bands. This is more follow-up 
telescopes than have been used for previous planetary microlensing discoveries. 

The light curve for this event had four distinct caustic crossing features, which were 
all observed with good-to-excellent sampling. The first is a weak caustic entry at HJD' ~ 
5006.05, which is observed by MOA. The second is a caustic exit at magnification A ~ 60 
at HJD' ~ 5006.6. This region of the light curve is covered by the CTIO, Danish, Liverpool 
and Wise telescopes. The next light curve feature is a strong caustic entry, which produced 
the light curve peak at A max ~ 205, at HJD' ~ 5006.96. The final caustic exit occurs shortly 
thereafter at HJD' ~ 5007.0 at a magnification of A ~ 180. This main peak covering the 
third and fourth caustic crossing has excellent coverage, observed by 16 telescopes. 

The images were reduced using several different photometry methods. The MOA data 
sets were reduced by the MOA Difference Image Analysis (DIA) pipeline (Bond et al. 2001). 
The /iFUN data sets except the CTIO H band and Bronberg were reduced by the MOA 
DIA pipeline and pySIS version 3.0 (Albrow et al. 2009), which is based on the numerical 
kernel method invented by Bramich (2008). The CTIO H band and Bronberg data sets were 
reduced using the OSU DoPHOT pipeline. The Danish data were reduced by the DIAPL 
image subtraction software (Wozniak 2000). The RoboNet and PLANET data sets were 
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reduced by pySIS version 3.0. The IRSF data set was reduced by the DoPHOT pipeline. 
The error bars for the data points are re-normalized so that x 2 P er degree of freedom for the 
best fit model is nearly one. 

All of these data sets are used for modeling except for the CTIO V and H band, the 
Canopus and SAAO / band, and the IRSF J,H,Ks bands. The CTIO V^-band, the Canopus 
and SAAO /-band, and IRSF J, H, i^g-band data sets do not have many observations and 
do not cover the planetary deviation region of the light curve. The CTIO //-band data 
was not used in the modeling because of a cyclic pattern caused by intrapixel sensitivity 
variations and image dithering. For our modeling of microlensing parallax effects, we have 
used a binned data set in order to speed up the modeling calculations. Note that we checked 
that an analysis with unbinned data gives the same results. 

3. Modeling 

Inspection of Figure 1 indicates that the event exhibits a number of caustic crossings, 
so we expect that this event, like most planetary microlensing events will exhibit significant 
finite source effects. The first step in modeling is therefore to measure the source color, which 
then enables us to determine the limb darkening parameters for the various light curves. 

3.1. Source Color 

Once a microlensing model is found, the dereddened source color and magnitude [I, (V — 
I)]o can be determined by comparing the instrumental values of these quantities to those 
of the red clump (Yoo et al. 2004). This is described in Section 4. However, before a good 
model can be found, the limb-darkening coefficients must be determined, which requires an 
estimate of the source color. This seemingly endless loop can be broken by making a model- 
independent measurement of the instrumental source color from a regression of V-band flux 
on J-band flux (and then comparing this value to the instrumental clump color). We find 
(V — I)o = 0.82, as reported in greater detail in Section 4. 
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3.2. Limb Darkening 

We adopt a two-parameter square-root law (Claret 2000) for the surface brightness of 
the source, 

S x (0) = S x (0) \l - c(l - cos 0) - d(l - Vcosti)] . (1) 

Here, c and d are the limb darkening coefficients, S\(0) is the central surface brightness of 
the source, and $ is the angle between the normal to the stellar surface and the line of sight, 
i.e., sin $=0/6*, where 9 is the angular distance from the center of the source. 

Based on the source color estimate of (V — I) = 0.82, the source is likely to have a 
G8 spectral type and an effective temperature of T e g- = 5475 K according to Bessell & Brett 
(1988). We use limb darkening parameters from Claret (2000) for a source star with effective 
temperature T cS = 5500 K, surface gravity logg = 4.5 and metallicity log[M/if] = 0.0 as 
presented in Table 1. 



3.3. Best Fit Model 

We search for the best fit binary lens model using a variation of the Markov Chain Monte 
Carlo (MCMC) algorithm (Verde et al. 2003) due to Doran & Mueller (2004) and Bennett 
(2010) that frequently changes the "jump function" in order to find the x 2 minimum more 
quickly. There are three lensing parameters in common with single lens events, the time of 
the closest approach to the center of mass to, the Einstein crossing time t^, and the minimum 
impact parameter uq. Binary lens models require four additional parameters: the planet-star 
mass ratio q, the binary lens separation d, which is projected onto the plane of the sky and 
normalized by the angular Einstein radius #e, the angle of the source trajectory relative to the 
binary lens axis a, and source radius relative to the Einstein radius p = 6*/6e- In addition, 
for each data set and pass band, there are two parameters to describe the unmagnified source 
and background fluxes in that band. 

We begin by conducting a very broad parameter search. The parameter search has 
been conducted by two independent codes. We perform 300 separate \ 2 minimizations with 
initial parameters distributed over the ranges —5 < logg < — 1, — 3 < logd < 0.3, in order 
to identify the parameter regimes of models that could explain the light curve. The initial 
parameters with logrf > 0.3 were not prepared because of the d <H- d~ l symmetry. We find 
that the only models consistent with the observed light curve have q ~ 10 -4 and d ~ 1, 
and that the best fit model has q = (3.95 ± 0.02) x 10~ 4 , d = 0.97537 ± 0.00007, and other 
parameters as listed in Table 2. The projected position of the planet is pretty close to the 
Einstein ring, and therefore d was well constrained. The light curves and caustic of this event 
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are shown in Figure 1 and Figure 2, respectively, which resemble Figure 8 in Wambsganss 
(1997). Here we assumed no orbital motion of the planet around the host star in our model. 
So the d and a are the average separation and angle during half a day when the source is 
crossing the caustics. The changes of these parameters due the orbital motion during this 
period could be same order or slightly larger than the nominal MCMC error of the average 
values given above. These changes does not affect the results in this analysis because they 
are much smaller than the uncertainty given in Section 7. 



4. Source Magnitude and Color 

The dereddened source magnitude and color can be estimated as follows. We locate the 
clump in the color magnitude diagram (CMD) of stars within 2' of the source star, shown 
in Figure 3, with the following method. The stars, which are J < 17 mag and (V — I) > 1.5 
mag, were used for the clump estimate. Among them, the stars within 0.3 mag from the 
clump centroid were picked up. Note that the clump in the first turn was assumed. Then, 
the mean magnitude of 7 and mean color (V — I) were calculated using the stars within 0.3 
mag and replaced as the new clump centroid. This was iterated until the clump centroid 
position is converged. Therefore, we find the clump as [I, (V — /)] c ium P = (15.31, 1.91). The 
best model source brightness and color are obtained as [I, (V — = (19.82, 1.69) from 
the fits. With a 0.05 mag correction due to blending by fainter stars in this crowded field 
(Bennett et al. 2010), this yields 

[J, (V - J)] s - [/, (V - I)] chunp = (4.51, -0.22). (2) 

We adopt the dereddened RCG magnitude M/ ,cium P = —0.25 and color (V — /)o, c ium P = 1-04 
from Bennett et al. (2008), which is based on Girardi & Salaris (2001) and Salaris & Girardi 
(2002). Rattenbury et al. (2007) find that the clump in this field lies 0.12 mag in the 
foreground of the Galactic center, which we take to be at Rq = 8.0 ± 0.3 kpc (Yelda et al. 
2010). Hence, the distance modulus of the clump is DM = 14.40. This yields a dereddened 
RCG centroid in this field of 

[/,(^-/)]cium P ,o = (14.15,1.04) . (3) 

Assuming that the source suffers the same extinction as the clump, we use the best fit source 
magnitude and color to obtain the dereddened values for the source, 



[I,{V-I)] s ,o = (14.15, 1.04) + (4.51, -0.22) 
= (18.66,0.82). 



(4) 



-12- 



A comparison of (V — I)s,o estimated by this method to 14 spectra of microlensed source stars 
at high magnification (Bensby et al. 2010) suggests that (V — I)s,o is determined with an 
uncertainty of 0.06 mag. For the uncertainty in ig , we estimate uncertainties of 0.08 from 
Rq, 0.05 from the Galactic bulge RCG centroid, and 0.05 from the Rattenbury et al. (2007) 
offset from the Galactic center, which when added in quadrature yields a total uncertainty 
of 0.11 mag. 

Equation (3) implies extinction of A\ = 1.16±0.11 and reddening E(V—I) = 0.87±0.08, 
which is consistent within the error with E(V — I) = 0.97±0.03 from the OGLE-II extinction 
map (Sumi 2004). 



5. Measurement of the Angular Einstein Radius, 6 E 

The sharp caustic crossing features in the MOA-2009-BLG-319 light curve resolve the 
finite angular size of the source star, and these finite source effects allow us to determine the 
angular Einstein radius #e and the lens- source relative proper motion fi re \ = #e/£e- Following 
Yoo et al. (2004), we use the dereddened color and magnitude of the source [I, (V— I)]s,o from 
Eq. (4). Next, we obtain the source angular radius using the source V and K magnitude. 
We estimate (V — K)q from (V — I)o and the Bessell & Brett (1988) color-color relations for 
dwarf stars, 

[K, (V - K)] Sfi = (17.67, 1.81) ± (0.14, 0.15). (5) 

We also estimate the K magnitude using IRSF data, it"s,o = 18.09 ± 0.42. This is consistent 
with but less accurate than the K magnitude estimated from (V — I)q. So, we use K 
magnitude estimated from (V — I)q. For main sequence stars, the relationship between 
color, brightness, and a star angular radius 8* was determined by Kervella et al. (2004) to 
be 

log(20*) = 0.0755(y -K) + 0.5170 - 0.2K, (6) 
which with K and (V — K) from Eq. (5) implies 

9* = 0.66 ± 0.06 /ias. (7) 

The fit parameter p = 6*/6e is source star angular radius in units of the angular Einstein 
radius. Thus, the angular Einstein radius 9e is 

Q 

e E = — = 0.34 ± 0.03 mas. (8) 
P 

Therefore, the source-lens relative proper motion \i is 

H = — = 7.52 ±0.65 mas yr^ 1 . (9) 
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6. Microlensing Parallax Effect 

The event time scale is not long, tg = 16.6 days, so one does not expect to detect the 
orbital microlensing parallax effect (Refsdal 1966; Gould 1992; Alcock et al. 1995). However, 
the very sharp third peak was observed simultaneously from Australia, New Zealand, and 
Hawaii, i.e., along two nearly perpendicular base lines of length, 0.36-R® and 1.25i? e , respec- 
tively. Therefore, there is some chance that these data will reveal a signal due to terrestrial 
microlensing parallax (Hardy & Walker 1995; Holz & Wald 1996; Gould et al. 2009). 

Microlensing parallax is usually described by the parallax parameter, 7Te, which is the 
amplitude of the two-dimensional microlens parallax vector, and the two components of this 
vector are denoted by 7Te.e and 7Te.n, which are the east and north components of the vector 
on the sky. The microlens parallax vector has the same direction as the lens-source proper 
motion, perpendicular to the line of sight. It is related to the lens-source relative parallax 
7r rc i and the angular Einstein radius #e (Gould 2000) by 

7TE = VT re l = 7T L - 7T S , (10) 

where 7Tl and 7Ts are the lens and the source parallaxes, respectively. 

Our initial search for microlensing parallax included both the orbital and terrestrial 
effect, as is necessary for a physically correct model. Our initial fits indicated a weak mi- 
crolensing parallax signal, so we searched for orbital parallax and terrestrial parallax signals 
separately, in order to determine which type of parallax signal is being seen and to test for 
possible systematic errors. We must also consider alternative model solutions due to the 
uq > -H- u < degeneracy first noted by Smith et al. (2003). As the model results listed 
in Table 2 indicate, orbital parallax can improve the fit x 2 by only A% 2 = 0.6, with two 
additional parameters, which is not, at all, significant. The best terrestrial parallax model, 
however, does give a formally significant x 2 improvement of A% 2 = 6.2, but this improvement 
decreases to Ax 2 = 6.1 for the best physical (terrestrial plus orbital) parallax model. With 
two additional parameters, this is formally significant at almost the 95% confidence level. 
Figure 4 shows the A% 2 contours for microlensing parallax fits to the MOA-2009-BLG-319 
light curve. 

The best fit parallax model has uo > and (7Te,e, ti"e,n) = (—0.15,0.15) ± (0.07,0.05), 
while the best fit u$ < model has a x 2 value that is larger than the best fit Uq > solution 
by 1.7 and only an improvement of A% 2 = 4.4 over the best fit non-parallax solution. Thus, 
the best no < model is neither a significant improvement over the best non-parallax model 
nor significantly worse than the best parallax model. We find that x 2 improvement for the 
best fit parallax model comes from Mt. John observatory (MOA-II 1.8m and Canterbury 
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0.6m) telescopes alone, with a total x 2 improvement Ax 2 = 7.3, while the contribution of all 
the other data sets is Ax 2 = —1.2 (i.e. the parallax model is disfavored). One would expect 
that x 2 should improve for the many other data sets, and the fact that it does not suggests 
that the parallax signal may not be real. 

If we assume that the scalar parallax measurement of 7Te is correct, then it implies that 
the lens system is located in the inner Galactic disk. Due to the flat rotation curve of the 
Galaxy, the stars at this location are rotating much faster than the typical line of sight to a 
Galactic bulge star. As a result, the direction of the parallax vector (which is parallel to the 
lens-source relative velocity) is most likely to be in the direction of Galactic rotation, which 
is ~ 30° East of North. This is similar to the direction of the parallax vector for the best 
ito < model, but it is roughly perpendicular to that for the u > model. So, the u > 
solution appears to be disfavored on a priori grounds. 

Because of the low significance of the microlensing parallax signal and the indications 
of possible systematic problems with the measurement of the parallax parameters, we will 
use only an upper limit on the microlensing parallax effect in our analysis. 



7. The Lens Properties 

We can place lower limits on the lens mass and distance with our measured angular 
Einstein radius, 9e, and our upper limit on the amplitude of the microlens parallax vector, 
7Te. The lens mass is given by 

Op 

M=— , (11) 

fv7T E 

where k = 4G/(c 2 AU) = 8.1439 mas Mq 1 . With our upper limit from the previous section, 
7Te < 0.5, gives a lower limit on the total mass of the lens system, M > 0.08M Q . This 
implies that the lens primary is more massive than a brown dwarf and must be a star 
or stellar remnant. From Eq. (10), this implies that the source-lens relative parallax is 
7r rc i < 0.1 7 mas. 

The vast majority of source stars for microlensing events seen towards the bulge are 
stars in the bulge, and the MOA-2009-BLG-319 source magnitude and colors are consistent 
with a bulge G-dwarf source. So, it is reasonable to assume that the source is a bulge star 
with a distance of D$ ~ 8.0 kpc. This implies that the lens parallax is 7Tl = vr rc i + 7Ts < 0.30 
mas, from Equation (10). The lens parallax is related to the distance by 7r L = 1A\J/D L , so 
a lower limit on the lens distance is > 3.33 kpc. 

An upper limit on the lens mass may be obtained if we assume that the planetary host 
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star is a main sequence star and not a stellar remnant. We can consider the blended flux seen 
at the same location of the source beyond the measured source flux from the microlensing 
models. If we attribute this blended flux to a single star, we can follow the reasoning of 
Section 4 in order estimate the dereddened magnitude of the blend star 

(/, V - /) b ,o = (17.78, 0.75) ± (0.12, 0.14) , (12) 

under the (conservative) assumption that the blend star lies behind all the foreground dust. 
We can now use this as an upper limit on the brightness of a main sequence lens star. From 
Schmidt-Kaler (1982) and Bessell & Brett (1988), we find an upper limit on the host star 
mass of M < 1.14M . 

As we found finite source effects in the light curve, we can break out one degeneracy of 
the lens star mass M, distance and velocity v. We calculated the probability distribution 
from Bayesian analysis by combining this equation and the measured values of #e and £e 
with the Galactic model (Han & Gould 2003) assuming the distance to the Galactic center 
is 8 kpc. We included the upper limit of microlens parallax amplitude. A constraint of 
the upper limit for blending light was also included for the lens mass upper limit. The 
probability distribution from a Bayesian analysis is shown in Figure 5. The host star is 
a K or M-dwarf star with a mass of Ml = 0.38lois ^© anc ^ distance Dl = 6.1^ 2 kpc, 
planetary mass M p = 50^24 M e and projected separation r± = 2. Ol2; 4 AU. The physical 
three-dimensional separation, a = 2.4lJ'g AU, was estimated by putting a planetary orbit at 
random inclination, eccentricity and phase (Gould & Loeb 1992). 



8. Discussion and Conclusion 

We have reported the discovery of a sub-Saturn mass planet in the light curve of mi- 
crolensing event, MOA-2009-BLG-319. This event was observed by 20 telescopes, the largest 
number of telescopes to participate in a microlensing planet discovery to date. The lens 
system has a mass ratio q = (3.95 ±0.02) x 10~ 4 and a separation d = 0.97537 ±0.00007 Ein- 
stein Radii. The lens-source relative proper motion was determined to be fi Te \ = 7.52 ± 0.65 
mas yr _1 from the measurement of finite source effects. A slightly better light curve fit can 
be obtained when the (terrestrial) microlensing parallax effect is included in the model, yield- 
ing an improvement of A% 2 = 6.1. This is very marginal statistical significance, and there 
are indications that systematic errors may influence the result. So, we use our microlensing 
parallax analysis to set an upper limit of 7Te < 0.5. 

The probability distribution estimated from a Bayesian analysis indicates that the lens 
host star mass is Ml = 0.38lo is ^0 w ith a sub- Sat urn-mass planet, M p = 50^4 M e and 
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the physical three-dimensional separation a = 2.4ljg AU. The distance of the lens star is 
.Dl = 6.1^2 kpc. The known microlensing exoplanets are summarized in Figure 6 and 
Table 3. MOA-2009-BLG-319Lb lies at ~ 2.3 times the distance of the snow line, which is 
estimated to be at a snow = 2.7 AU M/M Q . This is similar to the separation of other planets 
found by microlensing (Sumi et al. 2010). 

There is some indications of low level systematic deviations from the best fit model 
remaining in the light curve, near the third and fourth caustic crossing features (see the 
bottom panel residuals in Fig.l), which does not affect the results in this analysis. These 
systematic light curve deviations might be caused by orbital motion of the lens, a second 
planet, or systematic photometry errors. A more detailed analysis will be performed in 
the future when the adaptive optics images from the Keck telescope were reduced, and this 
analysis may shed more light on the mass and distance of the host star. 

The next few years are expected to see an increase in the rate of microlensing planet 
discoveries. The OGLE group has started the OGLE-IV survey with their new 1.4 deg 2 
CCD camera in March, 2010. This will allow OGLE to survey the bulge at a cadence almost 
as high as that of MOA-II, but with better seeing that should yield a substantial increase 
in the rate of microlensing planet discoveries. MOA also plans an upgrade to a ~ 10 deg 2 
MOA-III CCD camera in a few years, which will allow an even higher cadence Galactic bulge 
survey. The Korean Microlensing Telescope Network (KMTNet) is funded to dramatically 
increase the longitude coverage of microlensing survey telescopes. They plan three wide FOV 
telescopes to go in South Africa, Australia and South America. When these telescopes come 
online, we anticipate another dramatic increase in the microlensing planet discovery rate. 
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Fig. 1. — The light curve of planetary microlensing event MOA-2009-BLG-319. The top 
panel shows the data points and the best fit model light curve with finite source and limb 
darkening effects. The three lower panels show close-up views of the four caustic crossing light 
curve regions and the residuals from the best fit light curve. The photometric measurements 
from MOA, B&C, Auckland, Bronberg, CAO, CTIO, Farm Cove and LOAO are plotted as 
filled dots with colors indicated by the legend in the top panel. The other data sets are 
plotted with open circles. The data sets of /iFUN Bronberg and SSO have been averaged 
into 0.01 day bins, and the RoboNet FTN and FTS data sets are shown in 0.005 day bins, 
for clarity. 
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Fig. 2. — The caustic is plotted in solid curve for the MOA-2009-BLG-319 best fit model, and 
the dash line indicates the source star trajectory. The circle represents the source star size. 
The source star crosses the caustic curve four times, with peak magnification of A max ~ 205 
during the third caustic crossing at HJD' ~ 5006.96. 
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Fig. 3. — (V — I, I) color magnitude diagram of the stars within 2' of the MOA-2009-BLG- 
319 source using /iFUN CTIO data calibrated to OGLE-II. The filled triangle and square 
indicate the source and blend stars, respectively, assuming that the blended light comes from 
a single star. The filled circle indicates the center of the red clump giant distribution. 
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Fig. 4. — The contours of Ax 2 =l, 4, 9, 16 with orbital and terrestrial parallax parameters. 
The left panel is the result with uq > and the right panel is with uq < 0. The best fit 
result with uq > is better than uq < about A% 2 = 1.7. Furthermore, the best fit model 
with and without parallax is different with A% 2 = 6.1. 
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Fig. 5. — Probability distribution from a Bayesian analysis for the distance, Dj,, mass, Ml, 
and the physical three dimensional separation a. The vertical solid lines indicate the median 
values. The dark and light shaded regions indicate the 68% and 95% limits. The solid curve 
in the top panel indicates the mass-distance relation of the lens from the measurement of #e 
assuming D$ = 8 kpc. Note that D$ is not fixed in the actual Bayesian analysis. 
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Fig. 6. — Exoplanets as a function of mass vs. semi-major axis. The red circles with error 
bars indicate planets found by microlensing. The filled circles indicate planets with mass 
measurements, while open circles indicate Bayesian mass estimates. MOA-2009-BLG-319Lb 
is indicated by the gold-filled open circle. The black dots and blue squares indicate the 
planets discovered by radial velocities and transits, respectively. The magenta and green 
triangles indicate the planets detected via direct imaging and timing, respectively. The 
non-microlensing exoplanet data were taken from The Extrasolar Planets Encyclopaedia 
(http : / /exoplanet . en/). The planets in our solar system are indicated with initial letters. 
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Table 1. Limb darkening coefficients for the source star with effective temperature 
T cff =5500 K, surface gravity log g=4.5 and metallicity log[M/H]=0.0 (Claret 2000). 
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Table 2. The best fit model parameters with various effects, finite source, orbital and terrestrial parallax, and uq. 
The lines with V" list the la error of parameters given by MCMC. HJD'=HJD-2450000. Note that the u 
conventions are the same as Fig. 2 of Gould (2004). \ 2 value is the result of the fitting with 18 data sets, which have 
7210 data points. The model search with finite source and orbital parallax effects were done by a grid search. 
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0.07 





Table 3. Parameters of exoplanets discovered by microlensing. MOA-2007-BLG-400Lb has two solutions due to a 
strong close/wide model degeneracy, and details of the MOA-2008-BLG-310Lb parameters are discussed by Janczak 

et al. (2010) and Sumi et al. (2010). 





name 




Host Star Mass 


Distance 


Planet Mass 


Separation 


Mass estimated by 










M L (M ) 


D L (kpc) 


M P 


a(AU) 




OGLE- 


2003- 


BLG- 


235Lb 


n fi o +0.07 


r- o +0.6 

o.o _ a 7 


2.6 ±g:| Mj 


-0.8 


#e, lens brightness 


OGLE- 


2005- 


BLG- 


■071Lb 


0.46 ±0.04 


3.2 ±0.4 


3.8 ± 0.4 Afj 


3.6 ±0.2 


8e , tte , detection of the lens 


OGLE- 


2005- 


BLG- 


169Lb 


D AQ +0-23 
u -^ y -0.29 


O 7 +1-6 

-1.3 


13 ±| M ffi 


2 7 + 1 ' 7 

z -' -1.4 


#e, Bayesian 


OGLE- 


2005- 


BLG- 


390Lb 


22 +0 - 21 

u - zz -0.11 


6.6 ±i:8 


5.5 ±jjj M ffi 


2 6 + 1 ' 5 

z -° -0.6 


#e, Bayesian 


OGLE- 


2006- 


BLG- 


109Lb 

c 


n + 05 

u - 01 -0.04 


1.49 ±0.19 


231 ± 19 Me 
86 ± 7 M ffi 


2.3 ±0.5 

-1.0 


#E, tte 
#E, 7TE 


OGLE- 


2007- 


BLG- 


368Lb 


0-64 ™ 


5.9 ±H 


20 ±1 Me 


q q +1.4 
°-° -0.8 


#e, Bayesian 


MOA- 


2007- 


BLG- 


■192Lb 


0.084 tVoll 


0.70 ±°;? 2 


3.2 ±f:| M e 


0.66 ±°;^ 




MOA 


2007- 


BLG- 


400Lb 


0.30 ±g;S 


r o +0.6 

0,0 -0.8 


0.83 ±°;« Mj 


0.72 ±°;f 6 / 6.5 ±?J 


#e, Bayesian 


MOA- 


■2008- 


BLG- 


310Lb 


0.67 ±0.14 


> 6.0 


28 ±23 


1-4 i£S 


#e, Bayesian 


MOA- 


2009- 


BLG- 


319Lb 


0.38 t° fs 


o.i _ 12 


50 ±^ Me 


2-4 ± J 


#e, Bayesian 



